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Active optics Shack-Hartmann sensor: using spot sizes to measure 
the seeing at the focal plane of a telescope 
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ABSTRACT 

Real-time seeing estimation at the focus of a telescope is nowadays strongly emphasized 
as this knowledge virtually drives the dimensioning of adaptive optics systems and instru- 
ment operational aspects. In this context we study the interest of using active optics Shack- 
Hartmann (AOSH) sensor images to provide accurate estimate of the seeing. The AOSH prac- 
tically delivers long exposure spot PSFs - at the critical location of the telescope focus - 
being directly related to the atmospheric seeing in the line of sight. Although AOSH sensors 
are not specified to measure spot sizes but slopes, we show that accurate seeing estimation 
from AOSH images can be obtained with a dedicated algorithm. The sensitivity and compar- 
ison of two algorithms to various parameters is analyzed in a systematic way, demonstrating 
that efficient estimation of the seeing can be obtained by adequate means. 

Key words: Site testing - atmospheric effects - methods: numerical - methods: data analysis 



1 INTRODUCTION 

The evaluation of the seeing is paramount for selecting astronomi- 
cal sites and following their temporal evolution. Likewise, its esti- 
mation is fundamental for the dimensioning of adaptive optics sys- 
tems and their performance predications. This knowledge virtually 
drives instruments operational aspects at a telescope, and more em- 
phasis is made to develop and use accurate real-time seeing estima- 
tor at the focus of a telescope. 

The atmospheric seeing is commonl y measured by the dif- 
ferential image motion monitor (DIMM, ISarazin et al.J 1 19901) in 
most observatories, or by means of alt ernative seeing m onitors, 
e.g., Generalized Seeing Monitor (GSM. IZiad et al.feOOOh. Multi- 
Apert ure Scintillation Sensor (MASS, Kornilov and To kovininl 
1200 ll) . Being localized away from the telescope platform, such a 
device delivers seeing estimate that can significantly differ from the 
effective seeing seen at a telescope focii because of pointing orien- 
tation and/or height above the ground differences, or local seeing 
bias (dome contribution). The effect of the two latter will largely 
expand in the context of the next generation of telescopes: the ex- 
tremely large telescopes. 

On the other hand, the atmospheric seeing can be deduced at 
the focus of a large telescope from the width of the long-exposure 
point spread function (PSF) assuming successive corrections to ap- 
ply (turbulence outer scale, wavelength, and airmass). We remind 
the reader that the seeing is defined at A = 500nm for observations 
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at zenith. For this purpose, several flavors of images can be used at 
the critical location of the telescope focii: (1/) scientific instrument 
images, (2/) guide probe images, (3/) active optics Shack-Hartmann 
images. At the VLT, focal planes are equipped with an arm used for 
acquisition of a natural guide star. The light from this star is then 
split between a guide probe for an accurate tracking of the sky, and 
a Shack-Hartmann wavefront sensor used by the active optics to 
control the shape of the primary mirror. 

However seeing estimation from the full-width at half maxi- 
mum (FWHM) of a PSF strongly relies on the exposure time that 
must be long enough so that the turbulence has been averaged; en- 
suring that all representations of the wavefront spatial scales have 
passed through the pupil. This is thus dependent on pupil size 
and turbulence velocity, though it is commonly admitted that 30 
seconds average properly the turbulence, introducing significant 
FWHM biases otherwise. This already discards the use of the guide 
probe as exposure times are not longer than 50 milli-second, while 
scientific instrument images do not entirely comply with the real- 
time aspect of the seeing estimation. Instruments are affected by an 
observational bias: unavailability for a large range of seeing condi- 
tions, and are additionally affected by the telescope field stabiliza- 
tion. 

Active Optics Shack-Hartmann (AOSH) delivers continuously 
real-time images of long exposure spot PSFs (typically 45 seconds) 
at the same location of scientific instruments. AOSH images pro- 
vide simultaneously various data: slopes, intensities, and spot sizes. 
When short exposures are used, the information provided by both 
slopes and intensities (i.e., scintillation) can be used to retrieve Ctr 
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profile using correl ations of these data from two separated stars 
(Robe rtet al. J201 lh . In this paper, we propose to use the spot sizes 
in the sub-apertures to retrieve the atmospheric seeing in the line 
of sight. In other words, we propose to use the active optics AOSH 
sensor system of the telescope as a turbulence monitor to provide 
accurate seeing estimation directly at the telescope focus. For this 
purpose, careful assumption of the long-exposure PSF profile (to 
ensure precise FWHM evaluation), and accurate derivation of the 
seeing from the estimated FWHM are required. 

In this context, we investigate the performance and limitations 
of two different methods to estimate the seeing from long-exposure 
AOSH spot PSFs. The comparison and selection of the uppermost 
modus operandi is a byproduct of this work. Our study is carried 
out by means of extensive simulations, and using real data from 
the AOSH images obtained at the VLT. In Sect. 2 we describe two 
different algorithms developed to retrieve the seeing from long ex- 
posure images, and in Sect. 3 we present our simulation hypothe- 
sis, and discuss our results. In Sect. 4 real-data from the VLT AOSH 
database are re-analyzed and compared to synchronous image mea- 
surements obtained at t he VLT (UT4) with the instrument FORS2 
( Appenzell er. et aljl998l) . Finally in Sect. 5 we draw conclusions. 



2 EXTRACTING SEEING FROM SHACK-HARTMANN 
SPOT SIZES 

Accurate seeing estimation from the full-width at half maximum 
(FWHM) of a turbulence-limited long exposure PSF requires two 
conditions: (1/) careful assumption of the long-exposure PSF pro- 
file to ensure precise FWHM evaluation; (2/) accurate derivation of 
the seeing from the estimated FWHM. In the following, we treat 
these two aspects, and present in details the two dedicated algo- 
rithms. 

2.1 Long-exposure PSF profile 

The theoretical expression of a long-exposure PSF can be described 
through the expression of its optical transfer function (OTF) ob- 
tained by multiplying the telescope OTF, denoted T (f), by the at- 
mospheric OTF: 

r«CI) = exp[-0.5£VUf)], (1) 

where f is the angular spatial frequency, A is th e imaging wave- 
length, and DJ r) is the phase structure function dGoodmanll 1 9851 ; 
Rod dier]|l98ll) . This expression is miscellaneous and applicable to 
any turbulence spectrum and any telescope diameter. 
The analytic expression for the phase struc ture func t ion in the 
Kolmogorov-Obukhov model can be found in TatarskiJ dl96ll) and 
is expressed b y Dj,(r) = 6 .88(r/r ) 5/3 , where r is the Fried's co- 
herence radius dFriedl966l) . Finally, the long-exposure OTF can be 
expressed as: 

f(f) = r„(f) x exp[-3A4(At/r f 3 l (2) 

and long-exposure PSFs are accurately described by Eq. [2] In the 
case of a large ideal telescope with diameter D » r the diffraction 
term To can be neglected, while in the case of AOSH sub-apertures 
of size d it cannot (d ss r ). 

We note that Eq. [2] assumes nonrealistic behavior of the low- 
frequency of the turbulence model phase spectrum. It is firmly es- 
tablished that th e phase spectrum deviates from the pow er law at 
low frequencies dZiad et al. j200d : lTokovinin et al.]|2007t) . and this 
behavior is described in a first order by an additional parameter, 




Figure 1. Simulated AOSH image based on the VLT AOSH geometry. 

the outer scale Lq. The expression for the phase st ructure function 
(Dj,) with finite o uter scale L can be found in (e.g.. iTatarskiil 1 96lt 
Tokovinin 2002). As large-scale wavefronts are anything but sta- 
tionary, AOSH long-exposure spot PSFs can be in a first order de- 
scribed by Eq. [2] while a posteriori correction by L experimental 
estimate is mandatory. The examination of the Lq influence will be 
treated in the following. 

2.2 Seeing and FWHM 

The Kolmogorov-Obukhov model predicts dependence of the PSF 
FWHM £ on wavelength A and r : 

£ = 0.976 A/r Q . (3) 

Equation[3]is the definition of the seeing (assumed at A = 500 nm, 
and for observations at zenith). However, the physics of turbulence 
implies that the spatial power spectral density (PSD) of phase dis- 
tortions W$(y) deviates from the pure power law at low frequen- 
cies. A popular von Kar man (vK) turbulence model (e.g. . ITatarskiil 
ll96lUziadet"aL]|200Ch introduces an additional parameter, the 
outer scale Lq that describes the low-frequency behavior of the tur- 
bulence model phase spectrum: 

W,(V) = 0.0229 r V3 (|v| 2 + L 2 y n/6 . (4) 

Equation|4]is the definition of L , where v is the spatial frequency in 
m , The Kolmogorov-Obukhov model corresponds to Lq = oo. In 
the vK model, r describes the high-frequency asymptotic behavior 
of the spectrum. In this context, the outer scale of the turbulence 
(Lq) plays a significant role in the improvement of image quality 
(i.e., FWHM) at the focus of a telescope. The image quality is dif- 
ferent (and in some cases by a large factor, e.g., 30% to 40% in the 
near-infrared) from the atmospheric seeing that can be measured by 
dedicated seeing monitors, such as the DIMM. 
The dependence of atmospheric long exposure resolution on L is 
efficientl y predicted by a sim ple approximate formula (Eq.[5} intro- 
duced by Tokovinin (2002), and co nfirmed by means of extensive 
simulations (M artinez et al. J 1201 Pal o), where we emphasized that 
the effect of finite Lq is independent of the telescope diameter. The 
validity of Eq.[5]has been established in a L /r >20 and L /D^500 
domain, where our treat ment of the diffraction fa iled for small tele- 
scope diameters (D<lm, iMartinez et al.j|2010al) . 

FWHM ~ s Vl -2.183(r /L ) 0356 , (5) 
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As a consequence, to deduce atmospheric seeing e q (at 500nm) 
from the FWHM of a long-exposure PSF the correction implied 
by Eq.F5]is mandatory prior to airmass and wavelength correction. 
We will show that this does concern even the case of small size (d) 
AOSH sub-apertures, where Lo » d. 

2.3 Algorithms principle 

We consider two different algorithms to extract the FWHM of a 
long exposure AOSH image. These algorithms differ from their as- 
sumptions on the sub-aperture spot profile and the way sub-aperture 
diffraction is accounted for, while they rely on a common prelimi- 
nary step for the selection of spots retained for the analysis. 

2.3.1 AOSH spots selection 

For both algorithms, background estimation is performed on a cor- 
ner of the image without spots, and hot pixels are set to the back- 
ground. The cleanest and un-vignetted spots are selected in each 
frame for the analysis. These extracted spots are oversampled by 
a factor two, re-centered and averaged. The averaging reduces the 
influence of potential local CCD defects (e.g, bad pixels, bias struc- 
tures, etc.). In practice, the averaged spot is based on hundreds of 
selected spots. 

2.3.2 OTF '-based algorithm 

The first algorithm, hereafter Al, has been proposed by 
iTokovinin et al.J ( f2007l) and is based on the long-exposure spot PSF 
profile defined in Eq. [2] The modulus of the long-exposure opti- 
cal transfer function of the averaged spot is calculated and nor- 
malized. It is then divided by the square sub-aperture diffraction- 
limited transfer function To(f): 

7o(f) = (1 - [Aijd]) x (1 - [Af y /d]) (6) 

where d is the size of the AOSH sub-aperture. At this stage the cut 
of the T(f) along each axis can be extracted and fitted to the expo- 
nential part of Eq.[2]to derive a single parameter r , or equivalently 
Fourier transformed to derive the FWHM of the resulting spot PSF 
profile using a 2-dimensional elliptical Gaussian fit. The orienta- 
tion of the long and small axis of T(f) is found by fitting it with a 
2-dimensional elliptical Gaussian. 

2.3.3 PSF -based algorithm 

The second algori thm, hereafter A2, has been proposed by 
Noet he et al ., (2006), and is in use as a diagnostic tool at the VLT 
since May 2010. A2 relies on the assumption that AOSH spots can 
be described by the following rotationally symmetric PSF profile: 

F(u)ocF (u)® e xp[-(u/r ) 5/3 ] (7) 

where u is the spatial coordinate, F the diffraction limited PSF 
of the AOSH sub-aperture. The symbol ® denotes the convolution 
product. A2 relies on an approximate PSF profile based on the an- 
alytical expression of the OTF defined in Eq[2] The algorithm A2 
works as follow: an initial estimation of the FWHM (hereafter, 6) is 
derived by fitting the averaged spot profile to the exponential part 
of Eq|7] Then A2 approximately accounts for F a {u) by quadrat- 
ically subtracting the sub-aperture diffraction 8q = A/d from the 
estimated FWHM: 

FWHM « ^G 1 - d 2 . (8) 
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Figure 2. FWHM estimation from simulated AOSH images as function of 
the seeing obtained with Al and A2. 

This quadratic subtraction defined in Eq. [8] is justified by the as- 
sumption that the two terms of Eq.|7]can be approximated to Gaus- 
sian functions. 



3 CALIBRATION AND RESULTS 

In order to calibrate and compare these two algorithms, we test 
them with simulated AOSH images. To match our simulations 
to real situation, we assume identical conditions as the ones en- 
countered with the AOSH system of the VLT. Assumptions for 
the AOSH model and atmospheric turbulence are described bel- 
low. The wavelength considered through the study is 500nm except 
when the effect of the wavelength is analyzed. 

3.1 Simulation hypothesis 

3.1.1 Shack-Hartmann model 

Our simulations are based on a diffractive Shack-Hartmann model 
that reproduces the VLT AOSH geometry: 24 sub-apertures across 
the pupil diameter, 22 pixels per sub-aperture, 0.305" pixel scale, 
d = D/24 = 0.338 m. The validity of the AOSH model has been 
verified through several aspects such as the plate scale, spot sizes, 
slope measurements, and phase reconstruction. 

3.1.2 Atmospheric turbulence 

The atmospheric turbulence is simulated with 300 uncorrelated 
phase screens of dimension 3072 x 3072 pixels (i.e. 45 meters 
width). The principle of the generation of a phase screen is based on 
the Fourier approach: randomized white noise maps are colored in 
the Fourier space by the turbulence power spectral density (PSD) 
function, and the inverse Fourier transform of an outcome corre- 
spond to a phase screen realization. 

The validity of the atmospheric turbulence statistic has been 
carried out on the simulated phase screens: the value of the outer 
scale (Lo), Fried parameter (r ) and seeing of the phase screens have 
been confirmed by decomposition on the Zernike polynomials and 
variance measurements over the 300 uncorrelated phase screens. In 
addition the validity of the long-exposure AOSH image has been 
verified. In Fig.Q] we show an example of a simulated AOSH im- 
age. 
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Figure 3. Sub-aperture diffraction removal as function of the imaging wavelength applied on AOSH images generated at 500nm for Al (left) and A2 (right). 
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Figure 4. Al characterization: impact of the turbulence outer scale Lq (left) and the AOSH imaging wavelength (right). 



We have generated AOSH images through atmospheric turbu- 
lence with seeing conditions from 0.1 to 1.8" with 0.1" increment. 
Except when the turbulence outer scale is an open-parameter, L is 
always defined at 22 mete rs - VLT Paranal median value, see for 
instance lDaii. et al.l d2O10h . 



3.2 Results 

In this section we provide results from the comparison of Al and 
A2, while the best algorithm revealed will be further characterized 
against several important parameters. This calibration process is a 
necessary step prior to on-sky implementation of such numerical 
tool, and will provide insights for the re-analysis of the VLT AOSH 
database presented in Sect. 4. 



3.2.1 Algorithms comparison and selection 
- Response to the seeing: 

Figure [2] shows the FWHM extracted with both Al and A2 as 
a function of the seeing. The black dotted line corresponds to a 
pure equivalence between seeing and FWHM (i.e. Eq. [5] assuming 
L =oo), while the black dashed line corresponds to Eq. [5] with 
Lo=22m being in agreement with the statistic of the simulated 
atmospheric turbulence. The general trend of both algorithms is 



similar although A2 constantly provides smaller FWHM estima- 
tion than Al. In addition, Al provides a smoother response than 
A2 to the seeing; A2 exhibits irregularities in its behavior. 
Two regimes are observable: (1/) when the seeing is better than 
0.6" (Al) or 0.3" (A2) the FWHM is higher than the seeing, which 
is likely due to the under-sampling of the SH (pixel scale is 0.305"), 
(2/) when the seeing degrades further than 0.6" (Al) or 0.3" (A2) 
the FWHM is smaller and asymptotically converges towards Eq. 
[5] (L =22m). This last result is important as it demonstrates that 
although d <k L , the FWHM of a long-exposure PSF obtained 
with very small telescope diameters (here AOSH sub-apertures 
of ~30cm) does follow Eq. [5] We note that A2 convergence to 
Eq.[5](Lo=22m) is inaccurate for worst seeing conditions than 1.5". 



- Sub-aperture diffraction removal: 

Al and A2 differ from the way they account for the the sub-aperture 
diffraction. Al is based on the deconvolution of a square sub- 
aperture OTF (Eq. |6), while A2 makes use of a quadratic subtrac- 
tion of 6 (Eq. [SJ. To compare these two approaches, we assess 
the sensitivity of both algorithms on the sub-aperture diffraction. In 
practice, Al and A2 have been tested on the same set of AOSH im- 
ages simulated at 500nm, while the input imaging wavelength used 
to feed the algorithms for removing the diffraction A/d varies from 
400 to 800nm. Since the bandwidth center of the wavefront sensor 
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path can vary with the guide star type, pushing this test further with 
such a large range of imaging wavelengths (i.e., wider amount of 
diffraction) is here justified. 

Results are presented in Fig. [3] and show that Al and A2 
behave similarity for seeing conditions > 1.0", while for lower 
seeing by contrast to Al, A2 exhibits strong non reliable irreg- 
ularities. Subtracting quadratically the diffraction FWHM 6q as 
carried out by A2 is accurate as long as the effect of the diffraction 
is small relatively compared to the turbulence (d » r ), and 
therefore fails at large ro (good seeing conditions). In fact, the 
FWHM G varies from 0.24" to 0.49" when the wavelength varies 
from 400 to 800nm. The actual AOSH PSFs are a convolution of 
the atmospheric blur and diffraction, and since neither of these 
individuals broadening factors are Gaussian, calculation of the 
combined FWHM as a quadratic sum of individual contributions is 
not accurate. 

- Algorithm selection: 

Considering all the aspects treated previously, it appears that Al 
is more appropriate than A2. Therefore, in the following we will 
consider Al only for further characterization for the sake of clarity 
until Sect. 4, where real data obtained at the VLT with A2 will be 
re-analyzed. The next subsections concentrate on the dependency 
of Al to various parameters with the objective of allowing proper 
calibration of the algorithm. Nonetheless, it should be noted that A2 
demonstrates identical dependency to these parameters such that a 
calibration of Al or A2 can be carried out similarly. 



3.2.2 Turbulence outer scale 

In Fig. [4] (left) the extracted FWHM with Al as function of the 
seeing is shown for several outer scale (Lo) values: 10, 22m, and 
infinite outer scale (although likely about 200 meters due to the 
physical finite size of the phase screens used in simulation). It is 
observable that Al starts to be sensitive to the outer scale for seeing 
higher than 0.5". From results presented in Fig.|4](left) it is further 
established (see Sect. |3.2.TT > that AOSH sub-apertures are not small 
enough so that a spot FWHM can be approximated to eq, i.e., the 
AOSH spot FWHM measurement does depend on the outer scale 
and therefore follows Eq.[5] 



3.2.3 Imaging wavelength 

To analyze how the AOSH imaging wavelength impacts the estima- 
tion of the seeing, we generate AOSH images for various imaging 
wavelength covering the visible spectrum from 400 to 800nm. Re- 
sults are presented in Fig.|4](right) and indicate that accurate knowl- 
edge of the SH imaging wavelength is critical. For all wavelengths, 
Al response is constantly smooth and can therefore be efficiently 
calibrated. Al response always asymptotically converges towards 
Eq.[5](L =22m), where we remind the reader that r is wavelength 
dependent. The shorter the wavelength, the faster the convergence 
to Eq. \5j This behavior was already reported in iMartinez et al.] 
( 12010 a), where we examined the validity of Eq.|5]with wavelength 
in the context of telescope images. As noticed in Sect. 13.2.11 and 
generalized here for all wavelength, bellow 0.6" the impact of the 
AOSH pixel scale is observable. The reliability of seeing estima- 
tion fails at very good seeing conditions when seeing < 2x pixel 
scale (here 0.61"). 



3.2.4 Spot sampling 

To understand the asymptotical trend of the algorithm response to 
the seeing presented in the previous subsections, we analyze the ef- 
fect of the spot sampling on the FWHM measurement. Results are 
presented in Fig. [5] (left) where it is shown that for a given seeing, 
the estimation gets better when the sampling improves. FWHM at 
O'Vpixel sampling corresponds to theoretical values obtained with 
Eq.[5] A rough estimation indicates that 5 pixels per FWHM is re- 
quired for accurate measurement, though it can be calibrated and 
corrected in case of poor sampling. 

3.2.5 Field stabilization 

The telescope field stabilization removes the low frequency Tip-Tilt 
components generated by, among other things, wind shacking. As 
a consequence the field stabilization also removes the slow turbu- 
lence, which may reduce the FWHM of a PSF image. To assess the 
impact of the field stabilization we compare the FWHM measure- 
ments on simulated AOSH images where the Tip-Tilt contribution 
(in the full telescope pupil) has been completely removed (i.e., per- 
fect field stabilization with infinite bandwidth) to that of AOSH 
images left unmodified. This test has been carried out for several 
seeing conditions (0.1, 0.9, 1.8") and outer scale values (10m, 22m, 
and infinite case). No impact at all has been revealed but at a 0.01" 
level, which is negligible. The estimation of the seeing from the 
width of AOSH spots is therefore not sensitive to the field stabiliza- 
tion, which rises as a relevant advantage. The Tip-Tilt in the sub- 
aperture comes from the contribution of several low-order modes, 
which explains why the AOSH spot PSF FWHM is sensitive to the 
turbulence outer scale but not to field stabilization. 



3.2.6 Atmospheric dispersion 

Except for data taken at zenith, individual spot in the AOSH im- 
age can exhibit elongation in one direction due to uncorrected at- 
mospheric dispersion. Al allows the estimation of FWHM in two 
orientations: the small and the large axes. The algorithm has been 
tested on 200 real AOSH images obtained at the VLT-UT3 Nas- 
myth focus (not equipped with an atmospheric dispersion compen- 
sator, ADC) on May 12th 2010. Results of the FWHM extracted 
along the small and large axis are presented in Fig.|5](right) where 
the airmass is over-plotted. Results show that Al does differenti- 
ate the elongated axis, and the ratio between both axis follows the 
evolution of the airmass. The difference between the two axis can 
be as large as 0.3" which is substantial. The ability to distinguish 
small and large axis of the FWHM is therefore mandatory. We note 
that A2, by assuming a rotationally symmetric theoretical spot PSF 
profile (Eq.[7} does not allow to differentiate these two orientations. 

3.2.7 Detector PSF 

The diffusion of charges in the detector material before they are al- 
located to one pixel is observable as an artificial enlargement of the 
PSF. In most cases, the spatial response of the detector is not triv- 
ial to determine. Figure [6] (left) presents the impact of the detector 
PSF on the measured FWHM in the case of 1.0" seeing. We found 
that the detector PSF does enlarge the AOSH spot FWHM and that 
the effect can be significant; it starts to be substantial from 1 pixel, 
which is larger than what is usually encountered in scientific detec- 
tors, in particulars the ones at the VLT. 
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Figure 5. Left: Impact of the AOSH spot sampling. Right: Impact of the uncorrected atmospheric dispersion on the seeing estimation (real data). 
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Figure 6. Top: Impact of the detector PSF on the seeing estimation. Bottom: Impact of the signal-to-noise (RON=15 ADUs) on the seeing estimation. 



3.2.8 Signal-to-noise 

At the VLT the read-out-noise (RON) level of the AOSH is above 
15 ADUs, while typical signal is about 20000 ADUs (per sub- 
aperture). In figure[6](right) we present the measured FWHM in the 
case of 1.0" seeing as a function of the signal level (in ADUs) for a 
RON of 15 ADUs. It is shown that poor signal-to-noise (SNR) ratio 
image enlarges the FWHM which is evident for signal lower than 
10000 ADUs, while for signal higher than 10000 ADUs the FWHM 
measurements are roughly stable. The standard signal ADU level 
obtained at the VLT is therefore high enough to avoid any impact 
on the FWHM estimation, although it can be calibrated otherwise 
knowing the flux. 



4 APPLICATION TO THE VLT AOSH DATABASE 
4.1 Correction law and error budget 

To relate our simulations to real situation, we re-analyzed the VLT 
AOSH database obtained over the past year. At the VLT, FWHM 
estimates are recorded as diagnostic information since the commis- 
sioning of the telescope, while the A2 algorithm is used since May 
2010. As shown in the previous sections, A2 (similarly Al) requires 
to be calibrated. From Fig.[2]it is straightforward to fit the data of 
A2, and to derive a correction law to transform FWHM into seeing 



taking into account some of the critical parameters studied in the 
previous secti ons. Accounting f or the median outer scale value of 
Paranal (22m, iDali. et all201Ch . the under-sampling of the AOSH 
at the VLT, and the imaging wavelength, we found that the follow- 
ing equation provides accurate seeing estimate from A2 FWHM 
measurements: 



6 = 1.18 x (FWHM 184 - 0.15)" 2 



(9) 



From the sensitivity analysis of Al (likewise A2) to various 
parameters, we propose to draw here a classification of the critical 
parameters associated with an error budget. Based on the previ- 
ous systematic examination of all parameters impact on the seeing 
estimation, we can conclude that L , A, and the detector PSF are 
the major parameters to consider. All other parameters can be ex- 
cluded from the error budget, being either with negligible impact 
(field stabilization, SNR), with a known and constant value and thus 
calibrated (the spot sampling), or apart of the algorithm analysis 
(atmospheric dispersion). The Paranal median value for the outer 
scale is 22m. For the error budget of this term we consider the ± 
one sigma deviation around the median value 22m (11m and 42m). 
At the VLT, the bandwidth of the wavefront sensor path is centered 
around 550 nm in average. As it can vary with the guide star type, 
we consider a central wavelength of 550 nm ± 50nm. Hence we 
assume here a deviation between 500 and 600nm. For the detector 
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PSF term, since the value is unknown at the VLT, we consider that 
it is confined between our reference and ideal value and 1 pixel. 
We found that the error terms (cr) can be expressed as: 



ctl = 0.033 x eo, 
or., = -0.009 + 0.096 x e 



C del 



0.011 x eo 



(10) 
(11) 
(12) 



Assuming that all terms are independent, we can therefore write the 
error budget of Eq.[9]such as: 



'L 



+ cr,+ cr. 



(13) 



We note that for A2, by contrast with Al, an additional term should 
be considered: cr AD , standing for the atmospheric dispersion im- 
pact. From this budget error, we can estimate that all these inde- 
pendent error sources affect the seeing derived from AOSH images 
at ~10% level. Ideally, a drastic reduction of this error level is pos- 
sible through simultaneous and instantaneous measurement of Lq, 
and accurate knowledge of the guide star type. 



- F0RS2 


ond SH 


not corrected 


- F0RS2 


and SH 


corrected 




Figure 7. Simultaneous seeing estimation with FORS2 and the VLT AOSH 
database, while AOSH data are corrected by Eq.|9](blue dots), and left un- 
corrected (red dots). FORS2 data are corrected by Eq. f5](Lo). Best fit of 
each set of data is over-plotted (full red and blue lines), while the black 
full line represents the reference, i.e., the ideal perfect match of FORS2 and 
AOSH seeing estimates. 



4.2 Application to the VLT database 

We applied Eq.[9]retroactively on the VLT AOSH database obtained 
at UT4. In particular, we use a set of 6500 simultaneous FWHM 
measurements obtained since May 2010 with the AOSH (using 
A2) and the FORS2 instrument images. We remind the reader that 
FWHM estimation on long-exposure PSF from any scientific in- 
strument follows Eq.[5] We note that FORS2 data might be biased 
by several parameters, among of them the telescope field stabiliza- 
tion, and the accuracy of the FWHM extrac tion which is base d 
on the SExtractor software (Source Extractor, iBertin. et al.ll 19961) . 
Briefly, SExtractor is run on the reduced FORS2 images, and stars 
are identified based on Sextractor parameters, while their FWHMs 
are measured by the program through the use of a Gaussian fit ap- 
plied on the data. The FWHM is afterwards extracted from the 
estimated Gaussian profile. The examination of the reliability of 
FORS2 data is beyond the scope of this paper, and therefore the 
comparison between AOSH and FORS2 data must be here, care- 
fully considered. 

The correction to apply to the data is carried out as follow: 
(1) FORS2 FWHM data are corrected by Eq.[5] assuming Lo=22m 
(Paranal median value), while A2 FWHM are left uncorrected by 
EqlU (2) FORS2 FWHM data are corrected by Eq.|5](L =22m) and 
A2 FWHM data are corrected by Eq|9] Results are shown in Fig.|7j 
where the red dots correspond to (1) and the blue dots to (2). Fig- 
ure UJ demonstrates that FORS2 correction for L and AOSH cor- 
rection with Eq. [9] bring seeing correspondence in a better agree- 
ment. This is obviously expected as both FORS2 and AOSH de- 
liver FWHM estimate at the same location: the focus of the UT4 
telescope. This result further confirm the reliability of seeing ex- 
traction from AOSH data. 

The power-laws of both set of data are provided in the leg- 
end of Fig[7Jand are over-plotted for the sake of clarity. When the 
AOSH data are corrected by Eq. [9] it is observable that the slope 
of the power-law through the data is linear and equal to unity, by 
contrast to the situation where AOSH data are left uncorrected. 
Only a constant and homogeneous over-estimation of FORS2 see- 
ing remains. Such an overestimation reflects the presence of a bias, 
which likely finds its origin in the FWHM estimation carried out 
by SExtractor. For instance, long-exposure PSF profiles cannot be 



described by a Gaussian, and thus the calculation of the FWHM 
based on a Gaussian fit is not accurate. This typically induces an 
over-estimation of the seeing at ~ 10% level. 

We inform the reader that Eq.|9]can be used on the VLT AOSH 
database, being temporally valid from May 2010 up to now. From 
the commissioning of the VLT to May 2010, another correction law 
must be derived since another algorithm was in use. In addition, in 
the future, Al should replace A2 at the VLT, and thus Eq. [9] will 
require to be re-adapted accordingly. 



5 CONCLUSION 

Active Optics Shack-Hartmann (AOSH) sensor offers an advanta- 
geous possibility to estimate the seeing at the focus of a telescope: 
(If) it delivers long exposure PSFs, (2/) it is not affected by any 
observational bias (continuous real-time seeing estimation) by con- 
trast to scientific instruments, (3/) we show that it is not sensitive to 
the telescope field stabilization. 

The comparison of th e two proposed algorith ms gave advan- 
tage to the one proposed bv lTokovinin et all d2007l) which is based 
on the theoretical OTF expression of the long-exposure PSF. Al- 
though PSF-based and OTF-based algorithms exhibit similar gen- 
eral behavior, we found that Al provides a smoother response than 
A2 to the seeing and better account for the sub-aperture diffraction, 
and atmospheric dispersion. 

We show that the estimation of the seeing from AOSH im- 
ages is sensitive to several parameters but can be calibrated. We 
establish that even considering the small size of the AOSH sub- 
apertures (~30cm), it does depend on the turbulence outer scale 
L and therefore follows Eq. [5] We show that AOSH can be used 
to build statistics using median value of Lo, in a real-time fashion 
but relying on median Lo value, or requiring instantaneous mea- 
surement of Zfl, We note that in this context several independent 
campaigns h ave converged t o the median value of Lo of 22 m at 
Paranal (e.g jDali. etalfcoioh . 

The VLT AOSH database available since the commission- 
ing of the telescope does not provide seeing information but spot 
FWHMs. In this context, Eq.[9]is required to derive the seeing in- 
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formation. We remind the reader, that Eq.[9]is temporally valid until 
May 2010. With these considerations in mind, the re-analysis of the 
past years of the VLT AOSH database demonstrates a seeing cor- 
respondence in better agreement between AOSH and the FORS2 
imager. This result further confirms the reliability of seeing extrac- 
tion from AOSH images. 

Finally, the qualification and calibration of the algorithm Al 
are nearly completed and clear the path to its operational imple- 
mentation at the VLT for three months of test period starting in fall 
2011. Finally, we note that more emphasis is given to use closed- 
loop real-time AO instruments data in the near future to get the 
estimate of the seeing at the critical location of the telescope focus. 



ACKNOWLEDGMENTS 

The activity outlined in the paper was supported by the Euro- 
pean Commission, Seventh Framework Programme (FP7), Capaci- 
ties Specific Programme, Research Infrastructures; specifically the 
FP7, Preparing for the construction of the European Extremely 
Large Telescope Grant Agreement, Contract number INFRA-2007- 
2.2.1.28. We acknowledge Philippe Duhoux from ESO for his help 
with the VLT control software. 

REFERENCES 

Appenzeller et al., 1998, The Messenger, 94, 1-6 
Bertin et al., 1996, A&AS, 117, 393-404 



Dali et al., 2010, A&A, 524, A73+ 

Fried, D. L., 1966, Journal of the Optical Society of America, 56, 
1372+ 

Goodman, J. W., 1985, Statistical Optics by Joseph W. Goddman 
New York, NY: John Wiley and Sons, 1985 

Kornilov V. G. and Tokovinin, A. A., Measurement of the turbu- 
lence in the free atmosphere above Mt. Maoverset, Astronomy 
Reports 

Martinez, P., Kolb, J., Tokovinin, A., and Sarazin, M., 2010, A&A, 
516, A90+ 

Martinez, P., Kolb, J., Sarazin, M., and Tokovinin, A., 2010, The 
Messenger, 141, 5-8 

Noethe, L., Sarazin, M., 2006, ESO internal document 

Robert, C, Voyez, J., Vedrenne, N., and Mugnier, L., 2011, Pro- 
ceeding in Comprehensive characterization of astronomical sites 

Roddier, F, 1981, Progress in optics. Volume 19. Amsterdam, 
North-Holland Publishing Co., 1981, p. 281-376 

Sarazin, M., and Roddier, F, 1990, A&A, 227, 294-300 

Tokovinin, 2002, PASP, 114, 1156-1166 

Tokovinin, A., Sarazin, M., and Smette, A., 2007, MNRAS, 378, 
701-708 

Tatarskii, 1961, V.I. Wave propagation in a turbulent medium, 

Dover Publ., Inc., New York 
Ziad, A., Conan, R., Tokovinin, A., Martin, F, and Borgnino, J., 

2000, Applied Optics, 39, 5415-5425 



© 201 1 RAS, MNRAS 000,[j]-?? 



